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Abstract. - Competing interactions are often responsible for intriguing phase diagrams in cor- 
related electron systems. Here we analyze the competition of instantaneous short range Coulomb 
interaction U with the retarded electron-electron interaction induced by an electron-phonon cou- 
pling g as described by the Hubbard-Holstein model. The ground state phase diagram of this 
model in the limit of large dimensions at half filling is established. The study is based on dynam- 
ical mean field theory combined with the numerical renormalization group. Depending on U, g, 
and the phonon frequency tJo, the ground state is antiferromagnetically (AFM) or charge ordered 
(CO). We find quantum phase transitions from the AFM to CO state to occur when U — A ~ 0, 
where A characterizes the phonon induced effective attraction. The transition is continuous for 
small couplings and large phonon frequencies u>o and becomes discontinuous for large couplings 
and small values of too- We comment on the possible relevance of this work for Bai-^K^BiO;!. 
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Introduction. — A classical problem in condensed 
matter physics is that of the competing effects of the 
Coulomb repulsion of the electrons with the attraction 
generated by the electron-phonon coupling. This appears 
most prominently in the theory of conventional supercon- 
ductivity. There, the relatively weak attraction mediated 
by the phonons wins against the Coulomb repulsion with 
the help of retardation effects and the fact that the latter 
is renormalized to a reduced value at the phonon scale 
[T]. Only with this pseudo potential effect phonon in- 
duced superconductivity is credible in spite of the om- 
nipresent Coulomb repulsion. Beyond such weak coupling 
arguments the competition of instantaneous repulsion and 
retarded attraction has only started to be explored in 
recent years with the advent of reliable strong coupling 
methods. In many materials the low energy physics is 
dominated by different competing - often strong - interac- 
tions, which need to be studied simultaneously. For exam- 
ple consider strongly correlated systems such as the high 
T c cuprates [2], fullerides [3], manganites [J], and organic 
salts [3] and their intriguing phase diagrams. The purpose 
of this letter is to analyze the competition of instantaneous 
local Coulomb repulsion and the attraction mediated by 
phonons for different coupling strengths, and establish the 



resulting ground states, allowing for symmetry breaking in 
the magnetic and charge channel. 

As concrete model to study the competing effects we 
choose the combination the fundamental model of local 
electronic correlations, the Hubbard model with a 
classical model of electron-phonon coupling, the Holstein 
model [7j. Both models are probably oversimplified as 
to describe real materials in detail, but they can serve 
to obtain insights into dominant interaction effects. The 
combined Hubbard-Holstein (HH) Hamiltonian reads 



cj creates an electron at lattice site i with spin a, and 

b\ a phonon with oscillator frequency ljq, hi !<7 = <y ia Ci a . 
The electrons interact locally with strength U, and their 
density couples to an optical phonon mode with coupling 
constant g. Depending on U, g, u)q, and the filling factor 
n, the model is expected to display a variety of different 
phases including polaronic normal (N) state behavior, an- 
tiferromagnetic (AFM), charge (CO) and superconducting 
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Fig. 1: (Color online) The dynamic effective interaction U e g(u>) 
as a function of uj. 

(SC) order. It is important to map out these phases to see 
whether for a certain material an effective description in 
terms of the fundamental HH Hamiltonian is sensible. We 
focus on the case n = 1. 

Due to the large number of coupled degrees of freedom 
no exact solution of the HH model is available for the 
general case, and there are few analytical methods which 
respect the quantum nature of the phonons and allow for 
arbitrary coupling strengths U, g. For the pure electron- 
phonon problem perturbative schemes, such as Migdal- 
Eliashberg theory, can be very successful; strong Coulomb 
interactions can however not be treated. For the combined 
HH model there has been a lot of progress in recent years 
in one and infinite dimensions. In the d = 1 situation nu- 
merical methods can be applied with high accuracy and 
the phase diagram could be established [SHU]- To our 
knowledge for d > 1 general phase diagrams are still miss- 
ing. In the case of high dimensions the dynamical mean 
field theory (DMFT) [TJj becomes exact, and it can gen- 
erate non-perturbative solutions, such that it becomes the 
method of choice for our purpose of studying arbitrary 
coupling strengths. 

The competition of the interactions can be seen in con- 
cise form when integrating out the bosonic field, which 
yields an effective electronic interaction 



As depicted in Fig. [TJfor large uj the Coulomb repulsion U 
is dominant, ujq enters as a relevant energy scale at lower 
energy, and for |w| < loq the competition between the bare 
interactions is most important. 

The DMFT calculations deal with these competing inter- 
actions, and as a main result of this paper the ground state 
phase diagram of the infinite dimensional HH model at half 
filling emerges. This is likely to be good approximation for 
the three dimensional case, where collective excitations are 
not dominant. As shown in Fig. [3J the transition from an 
AFM state to a CO state occurs neaiQ {/ off = 0, where 

For the given parameters, there is a small tendency towards 




Fig. 2: (Color online) The phase diagram of AFM and CO 
in the [/-A-plane exemplary for cjq = 0.6t. The thin black 
line A ~ U gives a continuous transition and the thick line a 
discontinuous one. A nonzero CO/ AFM order parameter was 
found above/below the dashed line with diamonds/circles. The 
dashed lines with points gives the transition for phases with no 
long range order, a paramagnetic metallic (PM), bipolaronic 
(BP) and Mott insulator (MI). 

U e fi = U — A with A = 2g 2 /uj - For U c s < we find a 
CO and for U e g > an AFM ground state. The phase 
boundaries for the phases without long range order, the 
paramagnetic metal (PM), the Mott insulator (MI), and 
the bipolaron (BP) insulator, are also shown in Fig. [5] 
We see that only for large couplings the phase boundaries 
merge, but for smaller couplings other scales are impor- 
tant as analyzed earlier for ljo = 0.2 [13H15J . The decep- 
tively simple result for the phase boundary of the ground 
state phase diagram is expected in limiting cases, such as 
the antiadiabatic one, loq — > oo with A kept fixed, where 
as seen from ([2]) the HH model reduces to the Hubbard 
model with coupling constant U e g. For the general case it 
is not a priori clear that this applies, as this depends on 
the mutual renormalization effects of the couplings at low 
energy. However, our results indicate that this remains 
valid for general loq. In contrast to this universality the 
details and order of the transition depend on both, the 
couplings and the value of loq. 

The infinite dimensional HH model has received con- 
siderable attention and has been studied by DMFT. For 
instance, the formation of polarons |f 6] and the competi- 
tion between the Coulomb repulsion and the induced in- 
teraction [T7| were studied in the symmetric phase. Re- 
cently, the phase diagram of PM, BP and MI was estab- 
lished |I31I14| . where phases with long-range order were 
not allowed for. Allowing for AFM order, the effect of 
the Coulomb repulsion on the electron-phonon interaction 
was investigated [TB]. Here we extend these earlier cal- 
culations allowing for commensurate AFM and CO and 

A > U for the transition line, which can however hardly be resolved 
in the plot. 
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establish the full ground state phase diagram. In our cal- 
culations we also studied SC solutions, but we found that 
for finite ujo CO has lower energy at half filling, thus SC 
does not appear in the phase diagram. 

Formalism. — For our calculations we assume a bi- 
partite lattice with A and B sublattice, where the matrix 
Green's function can be written in the form 



i 



Ca.<t(w)Cb, ct (w) - e| V £ fc CaA u )J ' 

(3) 

with Ca,o-(w) = uj + fia.a- — E a . cr (a;), a = A,B, and k- 
independent self-energy [19] . For commensurate charge 
order, we have pa.u = P — h c , PB.a = p + h c and 
Sfl l( r(w) = Un - IU jCr (-a;)*, with n = (n A + n B )/2, 
n a = J2cr n a,cn where n a _ a = (n a>cr ). For AFM order, 
one has pa.o = P — Mb.ct = M + c^s, and the con- 
dition Eg i(7 (w) = EA-tr(w). We consider spontaneously 
ordered solutions where the symmetry breaking fields van- 
ish, h c , h s —J- 0. The matrix elements of the local Green's 
function G_(u>) can be calculated by integrating the matrix 
elements of © over the density of states, which we choose 
as semi-elliptic, po(s) = 2V D 2 — £ 2 /ttD 2 . We solve the 
effective impurity problem with the numerical renormal- 
ization group PHKH] (NRG) adapted to these cases with 
symmetry breaking. In these calculations we have chosen 
the discretization parameter A = 1.8 and we keep around 
1000 states at each iterations. The initial bosonic Hilbert 
space is restricted to a maximum of 50 states. This is jus- 
tified by the phonon occupation n p h which does not exceed 
values of n p h — 10 except when A is much larger than U. 

In the AFM case the A-sublattice magnetization, 
•^afm = in>A = (^A.f — nA,\)/2 serves as an order parame- 
ter. For CO we define $ co = {ua — l)/2. To identify the 
ground state of the system, we must compute the total 
ground state energy per lattice site, E tot = (H)/N, of the 
HH Hamiltonian (TT|) in the different phases. This gives 
generally, 

-Etot = Skin + Eu + E p h + Eg . (4) 

The first term is the kinetic energy, which reads 



(• r >) 



where pAB.k.o{u) = — ImG^B,fc,cr(w)/7r for the offdiago- 
nal Green's function in (j3|) and f(u>) is the Fermi func- 
tion. In the non-interacting case it can be evaluated an- 
alytically with pfc.(j(w) = S(u> — + p), and we find for 
half filling, n = 0, E° in = -4D/3tt, which for D = 2 is 
E® in ~ —0.8488. This can be used as reference energy. 
The interaction energies Ejj, E g can be calculated from 
expectation values. We have 



Ei 



Id 



(6) 



We distinguish between A- and B-sublattice values, which 
are equal in the AFM case, but not for the CO case. 
The third term for the total energy is generally given by 
the expectation value of the number of excited phonons, 
E p h = ujo^b). The detailed behavior of the various con- 
tributions to the energy in the different phases as well as 
static and dynamic response functions will be discussed 
elsewhere |22| . 

Results. - The half bandwidth D = 2t = W/2 is cho- 
sen as 2 in the following to set the energy scale. In Fig. [5] 
the phase diagram is shown for lu$ = 0.6 to exemplify the 
behavior neither too close to the adiabatic nor to the an- 
tiadiabatic regime. We carried out numerous calculations 
for other values of ujq and apart from the location of the 
point separating continuous and discontinuous transitions 
similar behavior was found for other choices. Notice that 
for physical optical phonons (u>q ~ 0.01VT — 0.2VF) this 
value is rather large and the later presented results for 
ujq = 0.2 can serve as a better guideline. Limiting cases 
of the phase diagram are known and easily understood 
on a qualitative level. Along the £/-axis, the pure repul- 
sive Hubbard model (at half filling on a bipartite lattice) 
is known to be AFM ordered at weak coupling [53], and 
this order is smoothly connected to the strong coupling 
Heisenberg AFM [241I25] . Along the A- axis, the pure Hol- 
stein model has a charge ordered ground state for g > 
and uo > |26) , where the limits of weak and strong cou- 
pling are smoothly connected. For finite U and g we find 
that the transition line is approximately given by the line 
A ~ U. For the cases of weak and intermediate coupling 
(U < W) the order parameters become very small in our 
calculation close to the line U e g = 0. An example is shown 
for fixed U = 2 in Fig. [3] (top) , where the order parameters 
are plotted as a function of A. 

Near A = U the ordering scale is very small (< 10 -3 ) 
and can not be resolved well in our DMFT-NRG calcula- 
tions, which then do not converge sufficiently well. The 
transition appears to be continuous, as all relevant re- 
sponse quantities change continuously. There are strong 
indications that it occurs directly from an ordered to an 
ordered state and no intermediate regime exists. We draw 
this conclusion by considering the behavior for different 
values of uio (see also Fig. 0]) . For n = 1 in the antia- 
diabatic case, luq — > oo, it is known that for any finite 
[/ cff > the system is AFM ordered [53] and for U cS < 
in the CO or SC state [57] . We have therefore a continuous 
transition from an ordered to an ordered state at U e g = 
on varying U or A. In the DMFT-NRG calculations we 
find that the smaller ujq is the larger the order parameters 
near the transition become (see Fig. 0]). Thus we conclude 
that this transition scenario persists for weak coupling and 
finite luq. Mean field calculations in the adiabatic limit, 
ujq — > 0, also support the picture of a direct transition 
between ordered states, however, the transition is discon- 
tinuous then. From NRG calculations in the normal state 
we can calculate a local effective quasiparticle interaction 
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Fig. 3: (Color online) The expectation values $afm and <E> CC , for 
U — 2 (top) and U — 5 (bottom) as a function of A. The total 
energy can be seen in the inset. 



U r by comparing the magnitude of the lowest two-particle 
excitation energy with twice the one-particle excitations 
energy at the fixed point [23]. U r changes sign approxi- 
mately at U e g = 0, which is consistent with a change of 
ground states there. To examine in detail what happens 
exactly at the transition at weak coupling and to analyze 
the critical properties requires an effective theory specifi- 
cally for this transition. 

For larger couplings we have included in Fig. [5] a dashed 
line (o) above which DMFT-NRG finds solutions with fi- 
nite $ co and another one (o) below which $ a fm is well 
finite. An example for this behavior is shown in Fig. [3J 
(bottom) for U = 5, where a much sharper behavior is 
seen. The calculation of the total ground state energy (see 
inset of Fig. [3J shows that also here the transition occurs 
approximately at U c ff — 0. A number of quantities such as 
the double occupancy (n-j-n^) show discontinuities at the 
transition. The total energy is a continuous function of A, 
but it displays a kink at the transition, such that a first 
derivatives will be discontinuous. The transition can be 
studied as a function of U for fixed A, and a very similar 
picture emerges for weak and strong coupling. For details 



we refer to \'22\ . 

As already alluded to, apart from the coupling strength 
the behavior near the transition also depends on the 
phonon frequency. In the adiabatic limit, loq — > 0, the 
phonons are represented by a static field, and one does 
not find a substantial modification of the ordered state as 
long as U e g does not change sign. When this happens the 
order changes abruptly and we find a discontinuous tran- 
sition. In the antiadiabatic limit we deal with an effective 




Fig. 4: (Color online) The expectation values <1> for fixed U = 2 
varying A. We show the results for ojo = 0.2 (full line), loq — 0.6 
(dashed line), and loq — > oo (dot-dashed line). 

Hubbard model with U e ff. The bare coupling becomes 
very small close to the transition, where it changes sign. 
This gives then behavior of a renormalized mean field the- 
ory with an exponential behavior of the order parameter, 
which goes to zero for U e g — > 0, and thus a continuous 
transition. Our calculations access the behavior between 
these limiting cases and we give examples for the behav- 
ior of $ for fixed U and varying A for ujq = 0.2,0.6 and 
luq —> oo in Fig. [5] We can see that the order for finite 
uiq is always larger than in the luq — > oo case. Our results 
give discontinuous transitions upon varying A for U > 4 
for uio = 0.6, whilst for ujq = 0.2 already for U > 3. 

More insight into the properties of the system are pro- 
vided by the electronic spectral function. Here we con- 
sider the (majority) local lattice Green's function Ga,^(u), 
which is given by the momentum sum of the diago- 
nal element of ((3]), and its spectral function pA^(uS) = 
— ImG,4,t( w )/ 7r - Pa,o{w) is symmetric at half filling in the 
normal state but becomes asymmetric in situations with 
symmetry breaking. Note that at half filling the spectra 
for minority spin in the AFM case, pA.i(w), and for the 
£>-lattice for the charge order, pb,\{oj), can be obtained 
from lo — > —lu. In Fig. [5] we show the spectral function for 
U = 2 and various values of A approaching the transition 
in the AFM and CO state for luq = 0.2 and w = 0.6. 
For comparison we have added the normal state spectral 
function for U = A as a dotted line in all cases. 
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Fig. 5: (Color online) The local A-lattice (majority) spectral 
functions pA,f(w) in comparison AFM (upper part) and CO 
state (lower part) for ujq = 0.2, 0.6 for fixed (7 = 2 and various 
A approaching the transition. The dotted line gives the spectral 
function in the normal state without symmetry breaking for 
U = A for comparison. 

In the AFM state for small A, the spectra fit well to the 
mean field description, where a square root divergence is 
found below the gap and square root increase above the 
gap. [2H The higher energy parts are little modified for the 
case of U = 2 apart from the broadening of the band edges, 
but no features which can be attributed to the phonons 
can be identified. When interpreting the spectra one has 
to take into account the broadening and the limited energy 
resolution of the NRG at higher energies, which limit the 
accuracy. On increasing A, the AFM order and magnitude 
of the spectral gap decreases. The electron phonon cou- 
pling is effective here in screening the repulsive [/-term. 
No polaronic features can be identified in the spectra as 
the coupling is fairly weak. For u>q = 0.6, the general form 
of the spectrum remains the same on increasing A, whereas 
for wo = 0.2 small features near the gap emerge. In this 
case one can see that the ordered solution and spectral 
gap persist for values of A up to ~ U, whilst for cj = 0.6 
the gap vanishes earlier on increasing A (see also Fig. 0] 
for $). 

In the CO state close to the transition, A > U, similar 
spectral characteristics of a weak coupling instability at 
the Fermi surface (w = 0) as in the AFM case can be seen 
both for loq = 0.2 and ujq = 0.6. In the limit loq — > oo 
we would find very similar behavior in the CO state on 
increasing A as in the AFM state on increasing U, but as 
illustrated around Eq. IJ5J) the phonon-mediated attraction 
has an w-dependent dynamic form, which leads to differ- 
ent behavior. One can identify a pronounced quasiparti- 
cle peak for interaction values near the transition which 
however becomes suppressed for larger values of A. This 
suppression can be partly due to the broadening in the 



NRG procedure as discussed in detail for superconducting 
solutions PHI- On increasing A a principal peak develops 
in the CO spectrum. Its position, is fairly well described 
by the spectral shift of fully polarized mean field theory, 
A m f = Un^_ a — 2A<I> co , for instance, A m f = 1.5 for A = 2.5 
and U — 2 (cf. CO spectrum for ojq = 0.2). 

We make a few general remarks concerning the phase 
diagram which has been established. It is remarkable that 
the phase diagram of the HH model in one dimension [3] 
is similar to Fig. [2] There one finds a Mott insulator 
with strong antiferromagnetic correlations, but no long 
range order, when U e s > 0, and a Peierls, charge ordered, 
insulator for U e g < 0. There is, however, a metallic re- 
gion with finite spin gap, but no charge gap in between 
these two phases [5l!9~ll 1 1 1 . For larger U this intermediate 
region shrinks until one finds a direct first order Mott- 
Peierls transition, similar to the present observation. A 
major difference with the high dimensional results is the 
real symmetry breaking in our case as well as the existence 
of the intermediate region, for which we find no indication 
here. 

We comment on the possible relevance of this work for 
the three dimensional compound Bai-^K^BiOa bearing 
in mind that there are still controversial issues. According 
to band theory the compound would be metallic for x = 
with a half filled Bi 6s band, but it shows CO [29] accom- 
panied by a lattice distortion, which hints towards a strong 
coupling of the electrons to the lattice. For x > 0.35, SC 
with fairly high T c « 30 K appears [515] . A strong coupling 
to an oxygen breathing mode at roughly 70meV is thought 
to play an important role (301131] . Thus, it has been pro- 
posed that the Holstein or HH model with effective band- 
width W ~ 4eV can provide an approximate description 
of the compound and account for the prominent features 
of the phase diagram including CO and SC |32H34| . Esti- 
mates for the dimensionless electron-phonon coupling con- 
stant, denoted by A here, based on density functional cal- 
culations vary substantially, A = 0.4 — 1.5 [31,35 37J. If 
we take A = Apo(0) — 1, such that A ~ 3.14i, then ac- 
cording to the phase diagram in Fig. [5]the experimentally 
observed CO state can be explained only if the residual lo- 
cal repulsion U is significantly less than A = 0.79W. This 
is consistent with the estimates for a residual positive U 
in an effective one-band model U/W ps 0.2 [35]. With the 
latter estimate for U even for A = 0.4 [37J a CO state is 
stable for x = 0. However, for such small A it is difficult to 
explain the high T c for superconductivity in the conven- 
tional theory. These issues, including appropriate values 
for A for different fillings, still have to be clarified by future 
research. Once suitable values are found, the approach 
presented here or refinements of it can make predictions 
for quantities such as the spectral gap or the transition 
temperature for the CO and SC state. 

A question raised in this work is why the transition be- 
tween the two ordered states occurs for U ~ A in such a 
generality. Both bare interactions produce opposite effects 
for the interactions of low energy quasiparticles. The fact 
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that the sign of the bare values U — A is reproduced for the 
effective quasiparticle interaction U r suggests that U and 
g 2 are renormalized in the same way. A reason for this 
is the specific form of the electron-phonon coupling to the 
local density and the phase space restriction at half filling. 
This is different away from half filling and for couplings to 
other modes such as of the Jahn- Teller type [3]. A renor- 
malization group study |10| has examined the competing 
interactions for the HH model in one dimension. Exten- 
sions to higher dimensions would be of great interest. 

Conclusions. — We analyzed competing interaction 
effects and established the ground state phase diagram of 
the infinite dimensional Hubbard-Holstein model at half 
filling with a transition line U ~ A. This implies that 
unlike superconductivity, where retardation is effective, 
phonon induced CO prevails only if A exceeds the Coulomb 
repulsion. However as seen in Fig. [3J a small excess can 
be sufficient to have a strongly ordered state. We also 
identified a point separating continuous and discontinu- 
ous transitions. A continuous quantum phase transition 
between different ordered state is rare and could be of in- 
terest for further study. We expect the presented phase 
diagram to be general as it holds for finite uq as well as in 
the adiabatic and antiadiabatic limiting cases. From the 
similarities with the one-dimensional phase diagram, and 
recent results in d = 2 [3H] (adiabatic limit), we conjecture 
that the form of the phase diagram possesses validity also 
in two and three dimensions. Extensions to finite tem- 
perature and other fillings would be desirable, where the 
competition can be studied in the case of superconductiv- 
ity. 

I wish to thank O. Gunnarsson, A.C. Hewson, G. San- 
giovanni, and R. Zeyher for helpful discussions, and W. 
Roller and D. Meyer for their earlier contributions to the 
development of the NRG programs. 

REFERENCES 

[1] Morel P. and Anderson P. W., Phys. Rev. , 125 
(1962) 1263. 

[2] Lanzara A., Bogdanov P. V., Zhou X. J., Keller 
S. A., Feng D. L., Lu E. D., Yoshida T., Eisaki H., 
Fujimori A., Kishio K., Shimoyama J. -I., Uchida T. 
N. S., Hussain Z. and Shen Z.-X., Nature , 412 (2001) 
510. 

[3] Gunnarson O., Alkali-Doped Fullerides: Narrow-Band 
Solids with Unusual Properties (World Scientific, Singa- 
pore) 2004. 

[4] Millis A. J., Nature , 392 (1998) 147. 

[5] Powell B. and McKenzie R. H., J. Phys.: Cond. Mat. 
, 18 (2006) 827. 

[6] Hubbard J., Proc. R. Soc. London, Ser. A , 276 (1963) 
238. 

[7] Holstein T., Ann. Phys. (N.Y.) , 8 (1959) 325. 
[8] Clay R. T. and Hardikar R. P., Phys. Rev. Lett. , 95 
(2005) 096401. 

[9] Hardikar R. P. and Clay R. T., Phys. Rev. B , 75 
(2007) 245103. 



[10 

[11 
[12 
[13 
[14 
[15 
[16 
[17 
[18 

[19 

[20 
[21 

[22 
[23 
[24 

[25 

[26 

[27 

[28 

[29 

[30 

[31 

[32 

[33 
[34 

[35 

[36 

[37 
[38 
[39 



Tam K.-M., Tsai S.-W., Campbell D. K. and Castro 
Neto A. H., Phys. Rev. B , 75 (2007) 161103. 
H. Fehske G. H. and Jeckelmann E., Europhys. Lett. 
, 84 (2008) 57001. 

Georges A., Kotliar G., Krauth W. and Rozen- 
berg M., Rev. Mod. Phys. , 68 (1996) 13. 
Jeon G. S., Park T.-H., Han J. H., Lee H. C. and 
Choi H.-Y., Phys. Rev. B , 70 (2004) 125114. 
Roller W., Meyer D., Ono Y. and Hewson A. C, 
Europhys. Lett. , 66 (2004) 559. 

Roller W., Meyer D. and Hewson A. C, Phys. Rev. 
B , 70 (2004) 155103. 

Roller W., Hewson A. C. and Edwards D. M., Phys. 
Rev. Lett. , 95 (2005) 256401. 

Sangiovanni C, Capone M., Castellani C. and 
Grilli M., Phys. Rev. Lett. , 94 (2005) 026401. 
Sangiovanni G., Gunnarsson O., Koch E., Castel- 
lani C. and Capone M., Phys. Rev. Lett. , 97 (2006) 
046404. 

Metzner W. and Vollhardt D., Phys. Rev. Lett. , 62 
(1989) 324. 

Wilson R., Rev. Mod. Phys. , 47 (1975) 773. 

Bulla R., Costi T. and Pruschke T., Rev. Mod. Phys. 

, 80 (2008) 395. 

Bauer J. and Hewson A. C, to be published (2010). 
van Dongen P. G. J., Phys. Rev. Lett. , 67 (1991) 757. 
Zitzler R., Pruschke T. and Bulla R., Eur. Phys. J. 
B , 27 (2002) 473. 

Bauer J. and Hewson A. C, Eur. Phys. J. B , 57 

(2007) 235. 

Freericks J. K., Jarrell M. and Scalapino D. J., 
Phys. Rev. B , 48 (1993) 6302. 

Micnas R., Ranninger J. and S.Robaszkiewicz, Rev. 
Mod. Phys. , 62 (1990) 113. 

Bauer J., Hewson A. C. and Dupuis N., Phys. Rev. 
B , 79 (2009) 214518. 

Hinks D. C, Dabrowski B., Richards D. R., Jor- 
GENSEN J. D., Pei S. and Zasadzinski J. F., Physica 
C: Superconductivity , 162-164 (1989) 1405 . 
Loong C.-R., Vashishta P., Kalia R. K., Degani 
M. H., Price D. L., Jorgensen J. D., Hinks D. G., 
Dabrowski B., Mitchell A. W., Richards D. R. and 
Zheng Y., Phys. Rev. Lett. , 62 (1989) 2628. 
Chainani A., Yokoya T., Riss T., Shin S., Nishio 
T. and Uwe H., Phys. Rev. B , 64 (2001) 180509. 
Vekic M., Noack R. M. and White S. R., Phys. Rev. 
B , 46 (1992) 271. 

Freericks J. R., Phys. Rev. B , 50 (1994) 403. 
Berger E., Valasek P. and von der Linden W., 
Phys. Rev. B , 52 (1995) 4806. 

Shirai M., Suzuki N. and Motizuki R., J. Phys.: 
Cond. Mat. , 2 (1990) 3553. 

Liechtenstein A. I., Mazin I. I., Rodriguez C. O., 
Jepsen O., Andersen O. K. and Methfessel M., 
Phys. Rev. B , 44 (1991) 5388. 

Meregalli V. and Savrasov S. Y., Phys. Rev. B , 57 
(1998) 14453. 

Vielsack G. and Weber W., Phys. Rev. B , 54 (1996) 
6614. 

Kumar S. and van den Brink J., Phys. Rev. B , 78 

(2008) 155123. 



p-6 



